SOLUTION OF THE GAMMA RAY BURST MYSTERY ? 

Nir J. Shaviv and Arnon Dar 

Department of Physics and Space Research Institute 
Israel Institute of Technology, Haifa 32000, Israel. 

■ Abstract 

Photoexcitation and ionization of partially ionized heavy atoms in 
highly relativistic flows by interstellar photons, followed by their ree- 
mission in radiative recombination and decay, boost star-light into 
beamed 7 rays along the flow direction. Repeated excitation/decay 
of highly relativistic baryonic ejecta from merger or accretion induced 
collapse of neutron stars in dense stellar regions (DSRs), like galactic 
cores, globular clusters and super star-clusters, can convert enough 
kinetic energy in such events in distant galaxies into cosmological 
gamma ray bursts (GRBs). The model predicts remarkably well all the 
CO '. main observed temporal and spectral properties of GRBs. Its success 

strongly suggests that GRBs are 7 ray tomography pictures of DSRs 
in galaxies at cosmological distances with unprecedented resolution: A 
time resolution of dt ~ 1 ms in a GRB can resolve stars at a Hubble 
distance which are separated by only D ~ 10 10 cm. This is equivalent 
to the resolving power of an optical telescope with a diameter larger 
q \ than one Astronomical Unit! 
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The origin of gamma ray bursts (GRBs), which have been discovered 35 years 
ago, is still a complete mystery [1]. Their observed isotropy in the sky, deficiency 
of faint bursts and the lack of concentration towards the Galactic center, in the 
Galactic disk and in the direction of M31, strongly suggest [2] that they are 
cosmological in origin [3]. A cosmological origin implies [4] that GRBs have 
enormous luminosities during short periods of time, 

L » lO 5O c4x0 6 Aft erg s'\ (1) 

where d = G^slO 28 cm is their luminosity distance, = 1O~ 6 6 ergs cm~ 2 s~ l is 
their measured energy flux, x is a bandwidth correction factor and AQ is the 
solid angle into which their emission is beamed. Moreover, the short durations of 
GRBs imply very compact sources. Relativistic beaming is required then, both 
in order to explain the absence of accompanying optical light and X-ray emission, 
and in order to avoid self opaqueness due to 77 — > e + e~ . All these considerations 
[4] seem to support the favorite cosmological model of GRBs; relativistic fireballs 
[5] formed by mergers of neutron stars (NS) or neutron stars and black holes 
(BH) in close binary systems [6] due to gravitational wave emission [7], or by 
accretion induced collapse (AIC) of NS and white dwarfs (WD) [6,8]. Indeed, the 
above considerations and the observed rate of GRBs favor NS mergers/AIC as the 
source of cosmological GRBs. Nevertheless, no mechanism has been convincingly 
shown to be able to convert a large enough fraction of their binding energy release 
into 7 rays and to explain simultaneously the complex light curves and spectral 
behavior of GRBs [4]. 

Most of the binding energy released in NS mergers/AIC is expected to be in 
the form of neutrinos, gravitational waves and kinetic energy of ejected material 
[9]. No mechanism is known which converts efficiently gravitational waves into 
7 rays. Neutrino annihilation [6] [yv — > e + e~), and neutrino pair production in 
strong magnetic fields (vy v — > ue + e~) near merging/collapsing NS cannot convert 
enough binding energy into e + e~7 fireballs or relativistic e ± beams which can 
produce cosmological GRBs, because of baryon contamination [10] and the low 
efficiency of these processes [11]. Thus, if NS-NS and NS-BH mergers, or AIC of 
WD and NS, produce GRBs, the production must proceed through conversion 
of baryonic kinetic energy of a highly relativistic flow (a fireball or a jet) into 
7 rays. Here we propose that repeated photoexcitation/ionization of the highly 
relativistic atoms of the flow by star-light in dense stellar regions followed by 
emission of decay/recombination photons which are beamed and boosted to 7 
ray energies in the observer frame (see Fig. 1), produce cosmological GRBs. We 
show that this simple mechanism, which has been overlooked, is able to convert 
enough baryonic kinetic energy of highly relativistic flows in dense stellar regions 
into cosmological GRBs. We also show that it predicts remarkably well the main 
observed properties of GRBs; their burst size, duration distribution, complex 
light curves and spectral evolution. 



2 



We assume that the ejected mass in NS mergers or in AIC of NS or WD 
contains high Z nuclei (e.g., iron-like) either from the crust of the neutron star or 
from nucleosynthesis of ejected neutron star matter when it expands and cools 
[12]. In view of the uncertainties in modeling merger/AIC of compact stellar 
objects [11], rather than relying on numerical simulations, we deduce the total 
relativistic kinetic energy release in such events from observations of type II 
supernova explosions, which are driven by gravitational core collapse into NS or 
BH. In type II supernova explosions, typically, 1OM are accelerated to a final 
velocity of v ~ 7000 km s -1 [9], i.e., to a total final momentum P ~ 5M e v. 
Since core collapse is not affected directly by the surrounding stellar envelope, 
we assume that a similar impulse, / Fdt ~ P, is imparted to the ejected mass 
in merger/AIC collapse of NS or WD. If the ejected mass is much smaller than 
a solar mass, AM <C M , then it is accelerated to a highly relativistic velocity. 
Its kinetic energy is given then by ~ Pc ~ 3 x 10 53 erg. About 1% of this 
kinetic energy must be converted into 7 rays in order to produce cosmological 
GRBs with ~3x 10 51 erg. 

The natural birth places of close binary systems are the very dense stellar re- 
gions (DSR) in galactic cores [13], collapsed cores of globular clusters (GC) and 
super star-clusters (SSC) [14]. These DSR are rich in very luminous stars and 
have very large photon column densities, iV 7 > L/ncRe . For instance, the core 
of our Milky Way (MW) galaxy has a surface brightness, S ~ 1.2 x 1O 7 L pc' 2 , 
typical photon energies, e ~ 1 eV, and a radius, R ~ 0.3 pc [14]. These val- 
ues yield, on average, iV 7 > 10 23 cm' 2 . A similar value is obtained for the 
nearby M31 galaxy. Larger values are possible for the unresolved super star- 
clusters. Much larger values are obtained for active galactic nuclei (AGN). The 
highly relativistic ejecta expands, cools and becomes transparent to its own ra- 
diation after a relatively short time [5]. Its partially ionized, highly relativistic 
atoms see the interstellar photons blue shifted to X-ray energies. These pho- 
tons photoexcite/photoionize the high Z atoms of the flow which subsequently 
decay/recombine radiatively by isotropic X-ray emission in the flow's rest frame. 
For a highly relativistic flow with a Lorentz factor r » 1, the isotropic X-rays 
are beamed along the flow, and as long as 9 < 1/T their energies, e x , are boosted 
to 7 ray energies, 

e 7 ^2re x /(l + r 2 £ 2 ), (2) 

where 8 is the angle of the flow direction relative to the observer. The differential 
fraction of the emission that is directed towards the observer is given by 

dI/d(cos9) w 2r 2 /(l + v 2 e 2 ) 2 . (3) 

The total baryonic kinetic energy which is converted into beamed 7 ray emis- 
sion by repeated excitation/decay can be estimated as follows: The total pho- 
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toabsorption cross section of an atom in an atomic state n into a state n' is 



7re 2 



>dv= f nn ,, (4) 

m e c 



where the integral is over the natural line width and where the oscillator strength 
/„„/ satisfies the sum rule, J2 n ' fnn' = Z, with Z being the number of electrons in 
the system [15]. If the typical energy of the interstellar photons is e ~ 1 eV, and 
if their energy is boosted in the flow's rest frame by T ~ 10 3 , then their typical 
photoabsorption cross section by Fe atoms is 



/ a mn ihdu a h hcZ 
fe 2m e c Ye 



3 x l(T iS cm 2 . (5) 



This cross section is larger by about seven orders of magnitude than the typical 
cross sections for production of 7 rays in the interstellar medium by inverse 
Compton scattering, bremstrhalung and synchrotron emission by electrons. Thus, 
if the ionized or excited atoms recombine/decay radiatively fast enough, then 
the typical energy of GRBs (burst size) from NS mergers/ AIC in GC is given 
approximately by 

£ 7 « aN^{Te/m A c 2 )E K « 3 x 10 51 erg, (6) 

where uia ~ 56m p . In its rest frame the relativistic flow expands against the 
external radiation field until its internal pressure equals the external pressure 
due to photoabsorption of the Lorentz boosted photons of the radiation field. 
Its relativistic atoms are ionized by the photoabsorption of the Lorentz boosted 
interstellar photons and than recombine radiatively and decay by line emission. 
The ionization/excitation rate adjusts itself to the recombination/decay rate. 
A detailed treatment [16] shows that the relativistic flow cools to rather a low 
temperature (eV), its Fe atoms are only partially ionized and it has a relatively 
large density which indeed result in fast enough radiative recombination and 
decay , except, perhaps, near AGN. However, if NS merger/AIC does occur near 
a bright AGN, where L ~ 10 12 L o and R ~ 10 17 cm, i.e., iV 7 ~ 10 29 cm~ 2 , 
then a large enough fraction of the kinetic energy of the relativistic flow can be 
converted into a GRB even by inverse Compton scattering [17,18]. Such a GRB 
will be extremely short and structureless. 

To calculate the detailed predictions of our model, we have constructed a nu- 
merical Monte Carlo code which simulates the formation of cosmological GRBs 
by highly relativistic flows in dense stellar regions [16]. The code employs the 
quantum mechanical cross sections for the relevant photo excitation/ionization 
and radiative recombination/decay processes. The core of the MW galaxy [13] 
was used for modeling the stellar environment (density of stars, stellar luminosi- 
ties and stellar temperatures) in a typical GC. Initial distributions of T which are 
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consistent with theoretical considerations (e.g., shock acceleration) and observa- 
tions have been used. Using the GRB simulation code we have found [16] that 
(a) the main properties of the simulated GRBs are not sensitive to fine details 
and (b) the calculated light curves and spectral behavior of simulated GRBs re- 
produce remarkably well those observed in GRBs [1, 19-21]. In particular, the 
simulated GRBs look undistinguishable from the observed GRBs. This is demon- 
strated in Fig. 2 which presents a simulated light curve of a GRB and its power 
spectrum, P(w) = | / L(t)exp(iwt)dt\ 2 , that has the universal power-law behavior 
P(w) ~ w~ 2 , which was found for real GRBs [18]. This simple universal behavior 
can be derived analytically [16,18] from our model. Also the other main observed 
properties of GRBs [1,19-21] follow from our model. Here we summarize briefly 
their approximate analytical derivations [16]. For the sake of simplicity we ne- 
glect here general relativistic effects (e.g., time dilation and energy redshift ) and 
assume that the explosions are spherical symmetric and occur at the center of the 
GCs, that the stars within a GC have a uniform spatial distribution, the same 
luminosity and the same effective surface temperature, and that the energy flux 
in the relativistic flow, E 2 driA/dE, is peaked around a Lorentz factor T. Then 
our model predicts that: 

1. GRB light curves are composed of a smooth background plus strong and 
weak pulses. Strong pulses are produced when the relativistic flow passes near 
luminous stars within the "beaming cone", i.e., stars at an angle 9* < 1/T relative 
to the line of sight from the explosion to the observer (see eq. 2). A multiple 
star system (binary, triplet, etc) produces a multipeak pulse with very short 
time-spacing between the peaks (spikes). Weak pulses are produced by boosting 
star-light from less luminous stars or by luminous stars outside the beaming 
cone. The smooth background is produced by boosting the background light in 
the beaming cone from all the other stars in the GC. 

2. The duration of a GRB reflects the spread in arrival times of gamma rays 
produced within the beaming cone (we neglect the short formation time of the 
relativistic flow because the dynamical time for NS merger/AIC is much shorter, 
typically few ms). It is given approximately by 



Hence, relativistic flows with T ~ 10 3 in MW-like GC (R ~ 0.3 pc) produce 
GRBs that last typically 30 seconds. Thus, approximately a one light-year path 
in a GC is contracted into 30 seconds 7 ray picture in the observer frame. 

3. The properties of the pulses from single stars depend on their locations, 
luminosities and spectra, and on the distribution of V in the flow. A pulse from 
a star at a distance and an angle 9* begins at a time ti ~ D*9 2 /2c after the 
beginning of the GRB (t = 0). Its duration (integrated intensity equals that due 
to the whole GC background) is given approximately by 



T ~ R/cT 2 . 



(7) 



T p ~ DM6* w b/cT ~ (47r/3) J RL,/7rcrL, 



(8) 
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where the integrated intensity from boosting photons from the star over impact 
parameters smaller than b equals that due to all the other stars in the GC. For 
main sequence stars with 1CT 1 < T*/T© < 10 3 , T = 10 3 and MW-like GC we 
obtain 10 ms < T p < 100 s. For typical early type stars with ~ 3 — 6L Q , 
T p ~ 0.3 — 0.6 s. Multiple stars yield multipeak pulses. For instance, the time 
difference between the two pulses from a binary star is given by 

T b « DM9*/c « d p /Tc, (9) 

where d p is the distance between the binary stars projected on the plane per- 
pendicular to the flow. Thus, T& < 500(c?A.u./r 3 ) ms, where c^a.u. is the binary 
separation in astronomical units. Since a large fraction of the stars are in close 
binaries, triplets, etc, a large fraction of the pulses have a multipeak structure. 

4. In a GRB from a GC with iV* stars uniformly distributed, the average 
number of stars within the beaming cone is n p ~ N*/4T 2 . Thus, for a MW-like 
GC and T ~ 10 3 , we expect, on average, n p ~ 3 stellar pulses. Their rate is 
dnp/dt ~ cA r *(r)/4r ~ 0.1 s -1 . It is independent of time if iV*(r), the number of 
stars within a radius r, is proportional to r, as seen in the GC of the MW [13]. 
Because n p is a small number, n p and the average time-spacing between strong 
pulses, AT, which is given by 

AT = T/rip ~ AR/cN* ~ 10s, (10) 

are predicted to fluctuate considerably around their average values. 

5. The total number of gamma rays in a stellar pulse arriving from the 
explosion is given by 

J F^dt « lOV 17 (L^L )/4^r^ lo m- 2 (11) 

where t* = D*/cT 2 = 10t w s is the peak intensity time of the stellar pulse and 
a = 10~ 17 (j 17 cm 2 . 

6. The durations of GRBs have a bimodal distribution which is a trivial con- 
sequence of their multipulse nature and the fact that, on average, T p AT inde- 
pendent of V: Multipulse GRBs have durations which are equal approximately to 
the sum of the time-spacing between their pulses. Consequently, T ~ SAT ^> T p . 
This produces a bimodal distribution which peaks around T p ~ 0.3 s for single 
pulse GRBs (see point 3) and around T ~ 30 s for multipulse GRBs. This is 
demonstrated in Fig. 3 for simulated explosions in a MW-like GC. 

7. Although the photoexcitation of partially ionized atoms produces line 
emission, the line emission is Doppler broadened into a continuum by the con- 
tinuous distribution of T of the atoms in the flow. Thus, the model predicts a 
continuous energy spectrum. The minimal Lorentz factor in the flow, r m , and 
the ionization state of the atoms cuts off emission in the observer frame below 
E m in ~ r m I ~ 5 — 25 keV, where I ~ few x 100 eV is the ionization potential of 
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the last bound electrons in the partially ionized Fe atoms [16]. It explains why 
GRBs are not accompanied by detectable X-ray or optical-light emission. 

8. If the spectrum of T in the flow has a power-law form, dnA/dT ~ T~ p for 
T > T m , then the energy spectrum of a strong pulse is given approximately by 

dn^/dE ~ cr(r)n e// (r)(rf J B/rfr)- 1 r- {p - 1 V(r^ + 1/2) (12) 

where a ~ r~( 1+<5 ) and n e ff is the average number of emitted photons per ex- 
citation/ionization. For T 3> 1/0*, one has [16] n e ff oc T, E ~ 2Fex, S < 1 
and then dn^/dE ~ £-(*>+ 5 ). For V < 1/^, one has [16] n eff « 1, £ ~ er 2 , 
5 ~ 0, and then dn-JdE ~ ^^(p +1 )/ 2 . Hence the spectrum of a pulse has a bro- 
ken power-law form. The break occurs when 2F0* ~ 1, i.e., around an energy 
E b ~ 2rex ~ ex/#* , and the power-index changes by ~ {p + l)/2 ~ 1.5. Since 
*i ~ D,9 2 j2c and « R, on average, £ b ~ 1/tf with 1/2 < a < 1. (For a GC 
with a uniform stellar distribution, the probability of -D* is proportional to D\ 
and consequently most of the stars have k, R.) 

9. The relative arrival times at a star of atoms with different V are given by 
t — ti D^cr -2 . Such atoms boost the star-light to an energy between e 7 ~ er 2 
and e 7 ~ 2exF. Therefore, the peak energy, E p = max E 2 (dn 1 /dE) ) decreases 
during a pulse approximately as ~ l/(t — U + <5tj) a i.e., 

E p ~ EiA*r Q (* - *i + ; 1/2 < a < 1, (13) 

where is proportional to the mean photon energy of the star-light and Sti is an 
added time broadening (a few ms) due to the time extension of the explosion and 
the extension of the region around the star where the boosting of the star-light 
takes place. E p is maximal right after the pulse begins and decreases monotonicly 
afterwards while the photon flux, dn^/dE, peaks at a later time, reflecting the 
1/r 2 behavior of the intensity of the radiation around a star and the distribution 
of T in the flow. Normal power-law spectra yield 7 ray fluxes during pulses that, 
on average, are time-asymmetric (fast rise and a slower decay) with longer pulses 
being more asymmetric. Pulses are narrower and their peak luminosities are 
shifted closer to the beginning of the pulse when viewed in higher energy bands 
(larger T) as demonstrated in Fig. 4. 

10. For a MW-like GC, iV 7 is not large enough to attenuate the relativistic 
flow. The kinetic energy of the the ejecta is injected into the galactic and in- 
tergalactic space. The general consequences of this injection will be discussed 
elsewhere. Here, we comment only on the high energy 7 ray production due to 
occasional collision of the ejecta with interstellar matter (stellar winds, planetary 
nebulae, etc) or with a molecular cloud with a sizable column density along the 
line of sight, in/near the GC (10% of the solid angle around AGN is covered by 
hydrogen clouds). Typical clouds, have R ~ 5 pc and M ~ 1O 4 M , yielding 
proton column densities of N p ~ 10 22 cm -2 . The total inelastic cross section of 
high energy iron nuclei on protons is a ~ 10~ 24 cm~ 2 . Thus, nuclear collisions 
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of the flow with a molecular cloud near the explosion will produce pions, and 
consequently, a burst of ~ 10 51 erg multi GeV 7 rays through tt° — > 27 decays 
(and neutrinos through — > fiu^ and /i — > ez/ e z/ M decays), with a power-law 
spectrum, dn^/dE ~ For p ~ 2, production of ~ 20 GeV 7 rays comes 

from 7r°'s produced mainly by nuclei with T ~ 200. Thus n° produced 7 rays of 
20 GeV are delayed, typically, by R/2cT 2 ~ 2/i, as was observed in the case of 
the 17 February 1994 GRB [22]. 

In conclusion, we have described a simple mechanism by which NS merg- 
ers/accretion induced collapse in dense galactic cores produce cosmological GRBs. 
The remarkable success of the model in reproducing all the main observed tem- 
poral and spectral properties of GRBs [16] strongly suggests that GRBs are 7 ray 
tomography pictures of DSR of distant galaxies with unprecedented resolution: 
A time resolution of dt ~ 1 ms in GRBs can resolve stars at a Hubble distance 
of 3000/i _1 Mpc ~ 10 28 /i _1 cm that are separated by only ~ cTdt ~ 3 x 10 10 cm. 
This is equivalent to the resolving power, ~ \/D, of an optical telescope with a 
diameter D > 2 x 10 13 cm, i.e. more than one Astronomical Unit! 

Finally, the proposed mechanism may play an important role also in other 
astrophysical gamma ray sources such as AGN, pulsars and other cosmic accel- 
erators. 
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Figure Captions 

Fig.l A schematic drawing illustrating the formation of a GRB by a highly rela- 

tivistic spherical flow in a dense stellar region. Most of the observed 7 rays are 

produced by the radiative decay of photoexcited atoms near stars within a cone 

with an opening angle 9 ~ 1/T along the direction to the observer. 

Fig. 2 A typical light curve of a simulated GRB and its temporal power spectrum. 

The straight line represents a w~ 2 power law behavior. 

Fig. 3 A duration distribution of simulated GRBs from a MW-like GC. 

Fig. 4 A simulated GRB from a MW-like GC viewed in different energy bands. 
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